function binomial()
%
% calculate and make sample data for a binomial release process.
%
%
% define parameters:
% m - quantal content
%
k=[0:15]; % range for binomial histogram
m=5; % quantal content
N=2500;  % number of observations




function [nk] = binomial(m, k, N)
nk = N*m.^k .* exp(-m)./factorial(k);





function [a] = alpha(tau, x)
alpha = talpha/ta * exp(-talpha ./ ta);
